Impact of van der Waals forces on the classical shuttle instability 
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The effects of including the van der Waals interaction in the modelling of the single electron shuttle 
have been investigated numerically. It is demonstrated that the relative strength of the vdW-forces 
and the elastic restoring forces determine the characteristics of the shuttle instability. In the case of 
weak elastic forces and low voltages the grain is trapped close to one lead, and this trapping can be 
overcome by Coulomb forces by applying a bias voltage V larger than a threshold voltage V u . This 
allows for grain motion leading to an increase in current by several orders of magnitude above the 
transition voltage Vu. Associated with the process is also hysteresis in the I-V characteristics. 

PACS numbers: 



I. INTRODUCTION 
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Nanoelectromechanical systems (NEMS) have seen a 
great upswing in interest over the past few years. NEMS 
show great promise for applications in areas such as 
charge detection, sensing, transistors and memory de- 
vices [E E S 0, 01 ■ In many NEMS devices the inter- 
play between electrical and mechanical degrees of free- 
dom is crucial for the transport of charge through the 
system. This is due to the Coulomb forces exerted 
by the charge carriers which may be strong enough to 
cause mechanical deformations of the systems. One sys- 
tem where such interplay drastically changes the charge 
transport properties is the single electron shuttle sys- 
tem |6j. There a metallic grain is embedded in an 
elastic material between two bulk leads. Since current 
through that system is accompanied by charging of the 
grain, interplay between the Coulomb forces and the 
mechanical degrees of freedom can lead to self oscilla- 
tions of the grain and electrons can be seen as shut- 
tled across the gap between the leads. Experimental 
shuttle related work has been reported on cantilevers 
that work as flexible tunneling contacts 0, , single- 



molecule transistors [ljj, colloidal particle systems 
and macroscopic shuttle systems [l^. Theoretical work 
has included different aspects of classical shuttle sys- 
tems d Il3l Il~il[l5l , some noise and accuracy consid- 
erations [lTJ 13 > different aspects of quantum mechan- 
ical shuttle models El El El El El El, and shut- 
tling of Cooper pairs in a superconducting shuttle sys- 
tem El Ell ■ F° r a recent review, see Ref. [27| . 

When down-scaling mechanical systems into the 
nanometer size range, the importance of short range sur- 
face forces, such as van der Waals (vdW) forces, adhesion 
forces and Casimir forces, increase causing problems with 
for instance stiction. To make realistic models of NEMS 
and to design NEMS-devices, or to make optimizations of 
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FIG. 1: Single electron shuttle, (a) A metallic grain of mass 
M and radius r placed between two leads separated by a 
distance L. The displacement of the grain from the center 
of the system is labelled X. The grain is connected to the 
leads via insulating elastic ligands. (b) Equivalent circuit of 
the system. The tunneling resistances and capacitances of the 
left and right junctions are Rl, Rr, Cl, and Cr. The leads 
are biased to the potentials Vl and Vr. 
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device designs, it is thus necessary to take such forces into 
account. In this work we report on a numerical investiga- 
tion of the effects of vdW-forces on the shuttle transport 
mechanism. The main result is that a new voltage scale 
comes into play, above which a grain trapped by the van 
der Waals force will be set loose due to the Coulomb 
interactions between the grain and the substrate. Asso- 
ciated with this process is a large increase in current as 
well as hysteretic behavior. 



2 



II. MODEL SYSTEM 

To be specific, we will here consider the impact of vdW- 
forces on a mechanically soft Coulomb blockade double 
junction. A schematic picture of this system is shown 
in Fig. A metallic grain of mass M and radius r 
is placed between two conductors (leads), separated by 
a distance L, via elastic, insulating materials (ligands). 
When a bias voltage V = (Vl — Vr) is applied between 
the leads, electron transport through the system occur 
by sequential, incoherent, quantum mechanical tunnel- 
ing of electrons between the leads and the grain. The 
system can, however, also lower its electrostatic energy 
by altering the grain position X. If the resulting electro- 
static force is comparable to the rigidity of the insulating- 
materials, grain displacement will occur. Thus, in this 
system, both electron transport and grain motion has to 
be taken into account. 



A. Electron transport 

In the notion of the equivalent circuit of Fig. ^> the 
electrical properties of the system are characterized by 
resistances and capacitances Rl, Rr, Cl, and Cr re- 
spectively. It is important that both the tunnel resis- 
tances as well as the capacitances depend on the grain 
position X. The resistances depend exponentially on X 
as Rl,r(X) — i?p R exp(±X/A) where A is the character- 
istic length scale for tunneling. The position dependence 
of the capacitances are given by 



Cl,r( x ) = 



a 



L,R 



1± 



where Cq ' R are the capacitances of the left and right 
junctions when the grain is at X = 0. The coefficients 
Al,r are typical capacitance length scales for the left and 
right tunnel junction respectively. 

Provided that the charging energy of the grain is larger 
than the scale of quantum fluctuations for all positions 
of the grain, only sequential tunneling need to be con- 
sidered and electrons will be localized in the leads or 
on the grain. This means Ec 3> fi/R(X)C(X), where 
Ec = e 2 /2C. Tunneling is then well described by the 
"orthodox" theory of Coulomb blockade |2^| from which 
it follows that the rates at which electrons tunnel through 
the left and right junctions are 
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Here (3 is the inverse temperature and AG L R (X,V,Q) 
is the decrease of free energy as the event (Q,Ql,r) — * 
(Q±e, Ql.rTs) occurs. The charges Ql,r and Q denote 
the charges accumulated on the left and right leads and 
the excess charge on the grain. The free energy decrease 



AG L R (X,V,Q) can be expressed using the equivalent 
circuit model of Fig. (2i| . 



B. Grain motion 

Considering the grain motion to be classical and 
one dimensional we have Newton's equation MX = 
FE X t.(X,V,Q) for the grain displacement. Whereas in 
Ref. Q the external force FE X t.(A, V, Q) acting on the 
grain was taken to be that of a simple damped har- 
monic oscillator driven by a term linear in Q, we are here 
considering a more realistic model including an elastic 
force Fk(X), a dissipative force, F 1 (X), an electric force 
F e (X, Q, V) and a vdW force F v dw{X) yielding the equa- 
tion of motion 

MX = F h {X) + Ey(X) + F e (X, Q, V) + F vdW (X). (2) 

The elastic force, Fk(X), has its origin in the compress- 
ible material between the grain and the leads. Since the 
grain displacement resulting from the electrostatic force 
may lie outside the range of validity of the harmonic ap- 
proximation, the nonlinearity of the elastic force has to 
be accounted for. In order to achieve this we have used 
a phenomenological potential that acts as a harmonic 
well with a usual spring force constant k for small dis- 
placements but that diverges as a Lennard-Jones poten- 
tial (12th power) at the positions Xz and Xr. The force 
is thus 



F k {X) 
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This elastic force is accompanied by a dissipative force, 



Fj(X) which is, as in 



compa 



phenomenologically mod- 



eled as a viscous damping term 

F y (X) = - 7 A. 

The electrostatic force, F e (X, Q, V), is given by the 
derivative of the electrostatic free energy G(X, V, Q) 



F e (X,V,Q) = -—G(X,V,Q) 

ldQr ^ldQ R 1 dV 

2^X VL+ 2^X VR -2 Q dX> 

and Q and Vl,r can be calculated from the equivalent 
circuit of Fig. ^p. 

The van der Waals force, F V( mr(X), between the grain 
and the leads can be derived from a Lennard-Jones in- 
teraction potential for the individual atoms j30| in the 
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FIG. 2: Elastic potential and vdW-potential plotted together 
with their sum. The whole graph is superimposed on a 
schematic figure of the system showing the grain and the 
leads. The dashed curve (red) shows the elastic potential 
arising from the the molecular links while the thin solid line 
(blue) is the vdW-force curve. The total potential is shown 
as a thick line (black). 



leads and the grain respectively. For a geometry with a 
metallic sphere close to a flat substrate one finds, 
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where H at t r ,rep are the Hamaker constants for the attrac- 
tive and repulsive parts of the force. This force has not 
been considered in this context before. 

Equation (J2J together with equation for the charge 
transfer rates describe the dynamics of the system. The 
complicated form of the potentials and forces in the 
present model makes a detailed analytical treatment 
cumbersome. However, there is no problem in treating 
the equations numerically and the ensuing results are eas- 
ily interpreted. Below we introduce a model system with 
a specific set of parameter values and investigate numer- 
ically the effect of incorporating adhesive surface forces 
on the system dynamics. The numerical investigation 
has been performed by doing Monte-Carlo style simu- 
lations incorporating a 4th-order Runge-Kutta routine 
when solving for the grain dynamics. 



III. RESULTS 

We have chosen numerical values for the model sys- 
tem with a typical state of the art experimental nano- 
electromechanical system in mind. We have considered 



a grain of mass M — 5 x 10~ 23 kg and radius r = 1 nm 
placed between two leads separated by a distance of 
L = 5 nm. The elastic potential was chosen to have a 
minimum in the center of the system, i.e. Xq — 0, 
and to diverge at a minimum grain-lead separation of 
0.5 nm which is just the statement that the ligands can 
never be compressed to less than a third of their size. 
The tunneling resistances and capacitances correspond- 
ing to zero grain displacement were R^ H = 10 9 tt and 
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F which corresponds to room temperature 



Coulomb blockade behavior. The tunneling length is 
material dependent but can be estimated to be of the 
order 1 A. The capacitance length scale is determined 
by the system size and geometry and is estimated to 
Al,r — 2.5 nm. For the vdW-interaction we have used 
the attractive Hamaker constant, H at t r = 4x 10~ 19 J [3(j 
and H rep was chosen from dimensional considerations, 
to be 10~ 72 Jm 6 . The focus on this work is the behav- 
ior of the shuttle instability as the relative strength of 
the elastic and van der Waals interactions change. We 
thus chose to consider an interval of the elasticity con- 
stant from k = 10~ 6 Nm -1 to k = 1 Nm -1 correspond- 
ing to small oscillation eigenfrequencies of the order of 
0.1 GHz to 100 GHz. This spans the interval of fre- 
quencies typical for NEMSpJ. The damping constant 
used was 7 = 10 -13 kg s _1 corresponding to quality fac- 
tors of approximately 0.5 x 10 3 through 0.5 x 10 6 which 
are also typical for those observed in NEMS 1]. In Fig. [2 
the vdW-potential and the elastic potential is plotted to- 
gether with their sum. The whole graph is superimposed 
on a schematic figure of the system showing the grain 
and the leads. The dashed curve (red) shows the elastic 
potential arising from the the molecular links while the 
thin solid line (blue) is the vdW-force curve. The total 
potential is shown as a thick line (black). It is clear that 
if, as is the case in our model, the insulating material 
does not allow the grain to come into direct contact with 
the leads the repulsive part of the vdW force does not 
come into play at all. It is also clear that the relation 
between the elastic force and the vdW-force determines 
the qualitative behavior of the total potential. Indeed, 
if the stiffness of the ligands is small the potential will 
have two stable minima whereas a stiffer configuration 
may have only one. The total potentials arising from 
different stiffnesses are shown in Fig. [3] 

In Fig. 0| we have plotted numerical results for the IV- 
characteristics as a function of the elasticity constant k. 
The figure is divided into three regions. In region I, the 
elasticity, k, is small and the potential has two stable 
minima; one in the vicinity of each lead. In region III, k 
is large enough for the potential to have a single min- 
imum in the center of the system. The potentials in 
Fig. |3| are representative of the different regions. The 
IV-characteristics in these regions differ markedly. Below 
we will elucidate this behavior with a mainly qualitative 
discussion of the dynamics, starting with the case when 
k is small, i.e. region I. Then moving on to discuss region 
III and finally we comment on the intermediate region, 
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FIG. 3: The total potential arising from elastic and vdW- 
interactions as a function of the grain position X shown for 
three different values of the ligand stiffness k. The solid 
line (black, k = 0.001 Nm" 1 ) correspond to region I in Fig. 01 
while the dashed line (red, k — 0.05 Nm -1 ) and the dash- 
dotted (blue, k = 0.3 Nm -1 ) correspond to the regions II and 
III respectively. 
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FIG. 5: (a) IV-characteristics for k — 0.001 Nm -1 at zero 
temperature. The inset shows a closeup of the transition re- 
gion. Here IV-curves for T = 0K for increasing (blue, thick 
solid) and decreasing (green, thick dashed) voltages are shown 
together with IV-curves for T = 300 K for increasing (black, 
thin dash-dotted) and decreasing (red, thin dotted) voltages, 
(b) Position of the grain as a function of time for V = 0.36 V 
in the case when the voltage is swept up. The grain sits asym- 
metrically in one of the potential minima caused by the strong 
van der Waals forces, (c) Position of the grain as a function 
of time for V = 0.36 V in the case when the voltage is swept 
down. Since the grain is already oscillating it does not get 
stuck in the potential minimum and the current is still large. 



FIG. 4: The current plotted as a function of the bias voltage 
V and the spring constant k. The figure shows that the tran- 
sition changes character from a very sharp transition in the 
region marked I to a more smooth transition in the region 
marked III. The intermediate region is marked II. In this 
graph all IV-curves correspond to increasing voltage sweeps 
which means that hysteresis effects are not shown here. 

labelled region II in Fig. 0] 

A. Small elasticity, region I 

In Fig. [SJl we have plotted the IV-characteristics for 
k = 0.001 Nm -1 . The most important characteristic is 
the big hysteresis loop in the IV-curve. The two lower 
panels (b and c) of Fig. [5] show the grain position as a 



function of time for low and high voltages respectively. 
For bias voltages around zero the the grain will settle in 
one of the two potential minima. Since the total cur- 
rent through the system is limited by the larger of the 
two tunnel resistances it will be strongly (exponentially) 
suppressed as compared to the situation when the grain 
occupies a more central position. As the voltage is in- 
creased the electrostatic force on the grain changes. The 
reason for this change is twofold; first, the electric field 
strength in the inter-electrode gap increases, and, second, 
increasing the bias will allow more charge to reside on the 
grain. When the voltage is further increased it eventually 
reaches a point V = V u where the grain can escape its 
locked position. In Fig. this happens at a bias voltage 
just above 0.4 V. Below this bias the grain cannot escape 
its locked asymmetric position, even though, as shown in 
Fig. [SJd, tunneling events can occur causing small move- 
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FIG. 6: Example of grain escape by a single pumping cycle. 
A grain initially carrying n = 4 extra electrons and at rest in 
the local minimum of the total potential cannot escape over 
the barrier at low temperatures. However, a charging and 
motion sequence as described by the arrows (1) through (4) 
provides the necessary energy to overcome the barrier. 



merits of the grain around the local potential minimum. 
Above V = V u , self oscillations of the grain appears, 
which causes an increase in the current by several orders 
of magnitude. 

The transition from the locked position can occur in 
several different ways. In the simplest case, in the ab- 
sence of thermal over-barrier activation and when the 
dissipation is large, the pumping due to grain motion 
inside the local minimum is very small, and the grain 
cannot escape until the minimum vanishes. However, if 
the minimum is shallow and the dissipation is not too 
high, correlation of charging and motion of the grain in 
the local potential minimum can lead to energy pump- 
ing and grain escape. An example of such a sequence 
of correlated events is illustrated in Fig. [fJl The figure 
illustrates a grain initially having n = 4 extra electrons 
and that is positioned in the local minimum of the total 
potential. The numbered arrows correspond to the fol- 
lowing events: (1) A tunneling event occurs changing the 
charge from n = 4 to n — 3. This change in charge is 
associated with a corresponding change in the potential. 
The grain is now not positioned in the minimum of the 
n = 3 potential. (2) The grain moves in the n = 3 poten- 
tial until a new tunneling event occurs. (3) A tunneling 
event changes the charge back to n = 4. (4) The grain, 
again in the potential for n = 4, now has enough energy 
to escape the potential well. 

The main IV-curve in Fig. [SJa, corresponds to the case 
of zero temperature. In this case, the maximum charge 
allowed on the grain, and hence, the most "shallow po- 
tential" experienced by the grain, is determined by the 
bias voltage due to charging effects (Coulomb blockade) . 
At finite temperatures, thermal fluctuations of the grain 
charge can cause the grain to become released at a lower 
voltage. In the inset of Fig. a close up of the transition 
region for both a finite (T = 300 K) and zero tempera- 
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FIG. 7: The IV-characteristics of the same system as de- 
scribed in Fig. but with different dissipation factors 7 to 
investigate features of the hysteresis. In (a) the dissipation 
is lowered to 7 = 10 -14 kg s _1 and we see that the hystere- 
sis is bigger than in Fig. |^| In (b) the dissipation is instead 
increased to 7 = 10~ 12 kg s" 1 and the hysteresis has disap- 
peared completely. 



ture is shown. As can be seen, the threshold were the 
grain escapes its locked position is considerably lowered 
for the high temperature case, while the remaining part 
of the curve remains unchanged. 

To understand the hysteresis appearing in the curve 
consider a voltage above V u . Provided that the dissipa- 
tion is not too high a sustained oscillation of the grain 
between the leads results. Whether or not this oscillation 
occurs depends, once the grain is released, on the relation 
between the work done by the electric field between the 
leads and the dissipation in the system 0, 0] . Lowering 
the voltage below V u once the grain is oscillating does 
then not imply that it will be stuck in one of the local 
minima of the potential. Hence, by changing the mag- 
nitude of the dissipation the size of the hysteresis loop 
varies. This is illustrated in Fig. [7| where a case with 
large hysteresis (7 = 10 -14 kg s _1 ) and a case with no 
hysteresis (7 = 10~ 12 kg s _1 ) are shown. It should be 
noted that the current in Fig. |3i and b differ by one or- 
der of magnitude. The dynamics of the latter case, when 
dissipation dominates the behavior, is mainly the same 
as in Ref. fl6| where it was further investigated. 
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B. Large elasticity, region III 

In the limit when the elasticity dominates the short 
range attractive forces only one minimum results from 
the elastic and vdW-potentials. This yields IV- 
characteristics belonging to region III in Fig. 01 The 
particular curve corresponding to k = 0.3 Nm" 1 is shown 
in Fig. [H^l. In this case, as the bias voltage is increased 
the shuttle instability develops from the one minimum in 
the potential. Sweeping the voltage down again reveals 
the absence of hysteresis. 

It should be noted here that, although the system op- 
erates in a regime analogous to the shuttle system in 
Ref. |(J, no steps are visible in the IV-curve even at zero 
temperature as would be the case if the elastic potential 
was purely harmonic. For a harmonic potential, increas- 
ing the bias leads to larger oscillation amplitudes at con- 
stant frequency, meaning that the grain comes closer to 
the lead. Since the tunnel resistance decreases exponen- 
tially with decreasing grain-lead distance the tunnel rate 
increases exponentially. Thus, increasing the amplitude 
at constant oscillation frequency implies that the grain 
charge has time to reach complete equilibrium with the 
lead. This in turn gives maximum step sharpness in the 
IV-characteristics |l7|. The presence of the strong an- 
harmonicity though, confining the grain to the center of 
the system, implies that with increasing bias the energy 
pumped into the vibrational degrees of freedom increases 
not only the amplitude but also the frequency. When the 
frequency increases the time spent by the grain in close 
contact with the lead decreases. Unless the amplitude 
is increased accordingly the grain charge does not have 
time to reach equilibrium with the lead, leading to the 
suppression of the steps seen in the IV-curves presented 
here. 

Although it is hard to tell from the curve in Fig. [S^ 
that shuttling is present, it is more obvious from Fig. [HJd 
and Fig. [SJ: where the position of the grain as a func- 
tion of time is shown for the voltages V = 0.2 V and 
V = 1.0 V. 



C. Intermediate elasticity, region II 

In the intermediate region the potential minima are 
not as deep as in region I (cf . Fig. [3J • An interesting 
phenomenon can occur in this regime; As the bias volt- 
age is raised the grain can sit in one minimum, then be 
pushed to the other, oscillate there for a short time, and 
then change side again. This behavior is shown in Fig.^ 
where the position of the grain as a function of time is 
plotted for k = 0.05 Nm' 1 and V = 0.1 V. 

Raising the bias, the grain switches side with increas- 
ing frequency (see Fig. [5J;) until it changes side twice in 
each oscillation period. Since the potential is flatter in 
this regime, the shuttle transition occurs at a lower volt- 
age than in region III. 
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FIG. 8: (a) IV-curve for k = 0.3 Nm~\ (b) Position of the 
grain as a function of time for V = 0.2 V. The grain oscillates 
in the single minimum in the potential. The position of the 
minimum changes with the number of charges on the grain 
due to the applied electric field, (c) Position of the grain as 
a function of time for V = 1.0 V. As the voltage is increased 
the grain starts to oscillate across the system. 



IV. CONCLUSIONS 

We have investigated the effects of van der Waals forces 
on the current transport through a shuttle system. Our 
results show that the relative strength between these 
forces changes how the transition from non-shuttle to 
shuttle transport occurs. In the region where the elas- 
tic force is small the transition is very sharp and there is 
a big hysteresis loop in the IV-curve, whereas in the re- 
gion where the elastic force is big the transition is softer 
and there is no hysteresis. Between the two regions is 
a crossover region where the transition voltage is signif- 
icantly lowered and where, for lower voltages, the grain 
alternatingly oscillates inside and between the minima of 
the potential. 
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FIG. 9: (a) IV-characteristics for k = 0.05 Nm" 1 . (b) Posi- 
tion of the grain as a function of time for V = 0.1 V. For low 
voltages the grain oscillates in one potential minimum until 
a tunnel event changes its charge. It is then is moved across 
the system to the other minimum where it oscillates until the 
charge changes again, (c) Position of the grain as a function 
of time for V = 0.4 V. For higher voltages the grain starts to 
oscillate across the system. 
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